Weak-Scale Hidden Sector and Energy Transport in Fireball Models of Gamma-Ray 

Bursts 
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The annihilation of pairs of very weakly interacting particles in the neibourghood of gamma-ray 
sources is introduced here as a plausible mechanism to overcome the baryon load problem. This way 
we can explain how these very high energy gamma-ray bursts can be powered at the onset of very 
energetic events like supernovae (collapsars) explosions or coalescences of binary neutron stars. Our 
approach uses the weak-scale hidden sector models in which the Higgs sector of the standard model 
is extended to include a gauge singlet that only interacts with the Higgs particle. These particles 
would be produced either during the implosion of the red supergiant star core or at the aftermath 
of a neutron star binary merger. The whole energetics and timescales of the relativistic blast wave, 
the fireball, are reproduced. 



The discovery of the afterglow jl) associated with some 
of the gamma-ray bursts (GRBs) and isotropy of the 
emitted radiation both support the view that the GRBs 
occur at cosmological distances with a redshift of the or- 
der of one. In this sense, GRBs offer a new observational 
tool for probing the early universe. 

There have been several suggestions concerning the 
generation mechanisms as well as the distributions of 
photons in the core of the astrophysical object, i. e., the 
central engine of the GRBs Among these, for exam- 
ple, the resonant production of gamma rays during the 
collision of two neutron stars is one possible mechanism 
||. Pertinent to a typical neutron star, the core region of 
the progenitor has a characteristic radius of Rq ~ 10km 
with roughly constant temperature fee 7b ~ 50MeV and 
matter density po ~ 10 14 g/cm 3 . The daily GRBs require 
a burst duration of 1 ms ^ r < 0.1 s with a total energy 
release of ~ 10 53 erg. 

Given the extensive parameters of the progenitor and 
the observed power spectrum, it is known that the 
gamma photons are trapped in the core region. Then 
one is to think about other possible agents to transfer the 
correct amount energy in the given time interval through 
the baryonic load. One possible alternative is neutrinos 
; however, the mixing of the flavour neutrinos with the 
sterile one is strongly suppressed in such matter densi- 
ties, and thus, the oscillation picture runs into difficulties 
[^| . Another alternative mechanism would come through 
the axions; however, the transferred power diminishes if 
the breaking scale of the Peccei-Quinn symmetry gets 
higher ||. In fact, the required axion mass cannot be 
reproduced in known axion models at all 

In this letter we work out a different scenario for trans- 
porting the energy outside the GRBs central engine. The 
basic agent of the process is a CP-even, presumably light, 
scalar particle, S, which has no baryonic charges. Such 
singlets have been proposed to take into account the 



nonobservation of the standard model (SM) Higgs parti- 
cle by increasing its invisible decay rate fel. For clarity 
of the discussion, we write down the effective Lagrangian 
describing the interactions between this SM-singlet and 
the photons as 



L int 9 \™ 2 S S 2 + X S S 4 



(1) 



where A^ is the photon, and m s and As designate the 
mass and the quartic coupling of the singlet, respectively. 
The original model || has an unbroken U(l) symmetry 
associated to the complex nature of this singlet. How- 
ever, for the purpose of this work the global phase of the 
singlet is not important so we take it real. In the frame- 
work of the hidden Higgs sector models the couplings 
above take the form 
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where X° s is the bare quartic coupling of the theory, and 
k mh stands for hSS coupling In writing these ex- 
pressions we neglected invariant masses in a given chan- 
nel compared to the Higgs boson mass. The vertex fac- 
tor A~ ( for /177 coupling is a slowly varying function of 
yfs, and it is of the order of one [jl0|. In all computa- 
tions below we will parametrize results in terms of A5 
and A 7 without going back to relations above. However, 
one keeps in mind that numerically A 7 ~ O(10~ 4 ) k. 

The effective Lagrangian (|j) describes a real scalar field 
interacting with itself and photons. The scattering pro- 
cesses following from this Lagrangian are shown in Fig.l. 
Here the diagram (a) represents 77 <-> SS scattering with 
which the conversion of the photons to singlets in the core 
of the GRBs burster occurs. Furthermore, conversion of 
singlets to photons outside (baryon-depleted region of) 
the GRBs central engine happens with the same process 
in backward direction. To have a description of the en- 
ergy transport through the strong baryonic load, it is 
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FIG. 1. Relevant scattering processes following from the interaction Lagrangian where the wavy lines represent the photon and 
dashed lines the singlet. The diagram (a) is responsible for conversion of photons to singlets and vice versa. On the other hand, the 
diagrams (b) and (c) determine the mean free path of the singlets. 



convenient to start with the conversion of the photons to 
the singlets. The relevant cross section reads as 

A 2 / 4m 2 \ 1/2 

<*n-«0-52;(i-^) (3) 

where s/s ~ lOOMeV is the total invariant mass of the 
annihilating photons. This process occurs in the core 
of the progenitor, and the produced singlet pairs move 
out through the surrounding baryonic loading. Before 
describing the journey of the singlets in the baryon load, 
it is convenient to compute the rate of energy conversion 
from photons to singlets. Denoting the four-momenta of 
the photons by fci = (u>i,k\) and k 2 — (w2,/c 2 ), the total 
amount of energy converted to singlet pairs per unit time 
per unit volume reads 

Q = J ^y^^y^ n ( Wl ) n ( W2 )( Wl + w 2KeiCr(77 -> SS) 

where n(uj) is the equilibrium Bose population of the 

photons and v re i = 1 — ki ■ fe/wia^ = s/(2u>iu>2) is the 
relative velocity of the two annihilating photons. After 
integrating Q over the volume of the core region, the total 
luminosity for photon-to-singlet conversion becomes 

C^ss « A 2 1 70 (Lj erg ■ sec" 1 (5) 



where 0(m 2 /(fcsT) 2 ) terms are neglected in comput- 
ing Q. This is a good approximation in the burster core 
where temperature is high enough. Assuming that this 
will be the luminosity observed on Earth, a compari- 
son with the GRBs standard candle luminosity requires 
A 7 ~ 10 -8 . As will be discussed later, A 7 is a loop- 
induced quantity in the weak-scale hidden Higgs sector 
models so that such a small number is naturally expected 
there. Despite all these, however, what is observed on 
Earth is £s,5^ 77 , that is, one has to convert the singlets 
back to photons to simulate the experimental conditions 
so that these naive bounds on A 7 may vary. 

The singlets, after being pair-produced by photon- 
photon annihilations in the core region, travel through 
the strong baryonic load towards the baryon-depleted 
region outside the GRBs burster. Since there is no in- 
teraction with the baryons they do not feel the baryonic 
load at all, and would move freely along radially out- 
ward trajectories were it not for their self-interactions. 
It is clear that the temperature of the host baryon dis- 
tribution does not effect the distribution and dynamics 
of the singlets as they can never come to thermal equi- 
librium with the baryons. In this sense, even if the star 
is quite cold with a small fraction of MeV temperature 
in the optically thin region, the singlets themselves can 
be quite energetic to generate the SS — * 77 reaction. 
Namely what singlets take out of the star is the energy 
accumulated in the gamma photons and this happens in- 
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dependent of the temperature and density distribution of 
the baryons. 

The singlet self-interactions are depicted in diagrams 
(b) and (c) of Fig.l. Both diagrams are generated by the 
singlet quartic coupling in the Lagrangian. A close in- 
spection of these diagrams reveal some important prop- 
erties. The SS — > SS scattering depicted in Fig. 1(b) 
preserves the number of singlets, is a contact interac- 
tion, and is kinematically operative for ^/s > 2m s . The 
SS — > SSSS scattering in Fig. 1(c), on the other hand, 
doubles the number of singlets, is a long-range inter- 
action with range ~ TnJ, and is kinematically allowed 
when yfs > Am s . The relevant cross sections are esti- 
mated to be 

a (Fig.l(fo)) = J£- » a (Fig.l(c)) ~ O (j^rA (6) 
1d7ts \ (2n)^s J 

so that cross section for SS — ► SSSS scattering is 
second order in X s and receives further phase space sup- 
pressions. 

The most important quantity describing the motion of 
the singlets towards the GRBs baryon-depleted region is 
their mean free path. Having no interactions with the 
baryons, the singlet mean free path would be infinitely 
long were it not for the singlet self-interactions depicted 
in diagrams (b) and (c) of Fig.l. The total mean free 
path obeys the relation 

where subscripts refer to the diagrams in Fig.l. There- 
fore, the total mean path over which singlets move freely 
is described by the larger of the individual contributions 
(for Fig. 1(b) and Fig. 1(c), respectively) 

i f d 3 p t , , / 4m2\ 1/2 

(b,c) = J (4^(^He* [l ~ ~f J O (8) 

where n(E t ) is the equilibrium Bosc population for the 
target singlet and v re i = s/ (2EE t ) is the relative velocity 
of the incident and target singlets with respective four- 
momenta p — (E,p) and pt = (E t ,p~t)- Here we take 
the phase space density of singlets as a Bose distribution 
ignoring the possibility of free streaming. In any case the 
resulting mean free path will be a conservative estimate 
of the actual one. 

As Eq.(||) suggests clearly, larger the cross section 
smaller the corresponding mean free path. Using the 
expressions for the cross sections in (^|) one can make the 
rough estimate 

„ „ ( E \ /50MeV\ /1(T 8 \ 2 
neglecting the terms O (ml/k^T^. 



This mean free path results solely from the self- 
interactions of the singlets, that is, it is the singlets them- 
selves which prevent their further flight. In particular, it 
is not the baryons that limit their motion so that it does 
not matter if £ happens to fall inside or outside the baryon 
load region; literally, singlets will stand still in a sphere of 
mean radius 0(£) measured from the core of the progeni- 
tor. To clarify this point further one recalls, for instance, 
the neutrino propagation. In that case the mean free 
path is determined by the neutrino-baryon interactions, 
and if it is outside the GRBs baryon load region neu- 
trinos get out of the astrophysical source otherwise they 
are trapped by the baryons and form a thermal neutrino 
sphere. Therefore, the formation of the singlet sphere 
follows only from their self-interactions. 

Due to their self-interactions in Fig. 1(b) and (c) sin- 
glets will form a thcrmalizcd cloud of particles whose 
number and energetics will change with a chain of such 
scatterings. As mentioned before, the 2^2 scattering 
in Fig. 1(b) preserves the number of singlets and plays an 
important role in restricting the singlets to have the finite 
path (|9j). The 2 — > 4 scattering in Fig. 1(c), however, is a 
long-range interaction and it modifies the number of sin- 
glets. Due especially to its long-range nature it is effec- 
tive everywhere in the singlet sphere, and causes largely 
separated singlet pairs to annihilate into four new sin- 
glets. This interaction, thus, increases the number of sin- 
glets and reduces the mean energy per capita. The result- 
ing cloud of singlets will thermalize themselves with their 
self-interactions with a temperature much lower than the 
burster core temperature. At any point inside the singlet 
sphere, there will be singlets coming from every direc- 
tion which is important in computing the energy accu- 
mulation in a given region. If singlets were moving along 
radially outward trajectories there would be a strong geo- 
metrical suppression factor for the energy deposition [jllj . 

As mentioned above, because of 2 — > 4 process in 
Fig. 1(c) total number of singlets increases, and thus, av- 
erage energy per singlet decreases. Similar to the elec- 
tromagnetic showers initiated by photons, a straightfor- 
ward computation of the total number of produced sin- 
glets as the mean energy per singlet drops from Eq to a 
critical energy E c via the 2 — > 4 scattering gives N — 
N [1 + (Eq — E c )/E c ] where N is the initial number of 
singlets given by the number of photons. In the problem 
at hand, Eq ~ fcslo and E c ~ ksT ~ 2m s where the 
latter follows from the kinematic blocking of 2 — > 4 scat- 
terings. Number of singlets per unit volume can be com- 
puted over a sphere of radius £: p = 3N/(4tt£ 3 ). bmce 
the singlet cloud eventually thermalizes it is convenient 
to use the usual Bose population for the number of sin- 
glets per unit volume per unit momentum state with an 
appropriate scaling to reproduce the singlet density p: 

j^n(E) (10) 
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where n(E) is the Bose population of the singlets at 
temperature T, and K is a constant for reproducing p. 

The next thing one needs for computing the electro- 
magnetic power generated by the singlets is the SS — > 77 
cross section 



a(SS -> 77) = 



\2 

A 7 



1 - 



Ami 



-1/2 



(11) 



to be compared with 77 — * SS cross section in (|^). 
Then computation of the luminosity proceeds in exact 
similarity with (|]) after replacing the cross section there 
by©: 



C 



SS->77 



i ) It) ^^ ss ' 



(12) 



where a factor of {T /T) 2 follows from (E /E c ) 2 . 
Css^-yy is computed over the volume of singlet sphere. 
Thus, the larger the mean free path the higher the lumi- 
nosity suppression factor. However, the higher the ini- 
tial (core) temperature the higher the luminosity emit- 
ted. Since both effects scale with the same power it is 
clear that no energy is lost in the global process, and 
thus the mechanism here suggested is able to effectibly 
carry out the overall energy released in the GRBs cen- 
tral engine as in the context of fireball models. Once the 
energy is out of the baryon loading, a relativistic blast 
wave of electron-positron pairs and radiation is formed, i. 
e., a fireball, which cleans away the burster environment 
pushing out a rather matter free debris to produce the 
observed burst when colliding with and external interstel- 
lar medium. The no-time-delay feature of this scenario 
makes it suitable to explain bursts with both rapid rise- 
time and prompt afterglows as observed in GRB990123. 
The potential role of this mechanism in triggering GRBs 
from supermassive star explosions is an issue currently 
pursued, jlj] 

To conclude, we have investigated the viability of the 
hidden Higgs sector models of turning realizable the en- 
ergy transport in GRBs in the context of fireball scenarios 
As the discussions in the text show, in such mod- 
els transport of the energy from the core to outside is 
quite efficient, and the resulting luminosity agrees with 
the astronomical observations. The particle physics sce- 
narios with neutrinos and axions are not as efficient as the 
present model due to the suppresion of the neutrino mix- 
ing angle and smallness of the axions mass, respectively. 
In the present model, conversion processes are mediated 
by the Higgs particle. On the other hand, transport of 
the energy from the core to outside is done by the singlets 
having a rather large mean free path compared to neutri- 
nos. In the present scenario singlets are light enough to 
be pair-produced by the photon annihilations, and such 
a light singlet does not contradict with the present col- 
lider data as it affects the precision observables at two 
and higher loop levels S. 



The model employed here has two free parameters, the 
coupling constants A 7 and As, both hidding the explicit 
dependence on the singlet mass m s and Higgs mass mh- 
The first parameter is constrained to be around A 7 < 
1CP 8 by the current GRBs BATSE observations, and its 
value will be measured at the 77 mode of the TESLA 
collider p5| . As a final remark, the potential realization 
of this scenario in the astrophysical sources triggering 
GRBs would render it a viable pathway for testing some 
of the extensions of the standard model of particle physics 
introduced to account for the overall GRBs observational 
properties: energetics, timescales, spectra, etc., such as 
the one being suggested here. 
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